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ABSTRACT 


Fibrous composites are fast replacing conventional 
engineering materials in various industries. Understanding the 
mechanism of fracture for these orthotropic materials is a 
problem. The theory of orthotropic photoelasticity can be 
employed for an efficient analysis of fracture phenomena in such 
materials. However, usage of this tool of experimental mechanics 
is possible only if a suitable model material is available. 

It is in this context that a model material exhibiting 
orthotropic properties, for photoelastic studies has been 
developed. Its optical propei'ties have been found to be much 
better than those of the materials developed earlier. 

Further, the principles of orthtropic fracture mechanics 
have been employed to determine the value of stress intensity 
factor for unidirectional laminates. Knowledge of this parameter 
is very important, since it can be used as a tool to understand 
the mechanism of fracture in composites. The values of S.I.F for 
the unidirectional laminate were found to be reasonably correct. 

It is expected that the newly developed model material and 
principles of orthotropic photoelasticity will meet the growing 
demand for an efficient stress and fracture analysis of fibrous 


composites. 
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1 INTSODDCTION 

1.1 Importance of Composite Materials: 

The word “composite" means 'consisting of two or more 
distinct parts’ [1]. Thus a material having two or more distinct 
constituent materials may be considered a composite material. It 
is only when the constituent phases have significantly different 
physical properties and thus the composite properties are 
noticeably different from the constituent properties that we have 
come to recognize these materials as composites. 

Composite materials have slowly emerged as easily workable 
engineering materials and are now quite commonly available and 
used, globally. Modern engineering industry finds application of 
these materials in areas as different as low technology 
applications like panels for doors to high-tech applications viz. 
space programs and aircraft industry. 

The success of fiber composites results from the ability to 
make use of the outstanding strength, stiffness, and low specific 
gravity of fibers such as glass, graphite, or Kevlar. When 
outstanding mechanical properties are combined with the unique 
flexibility in design capability and ease of fabrication that 
consposites offer, it is no wonder that their growth rate has far 
surpassed other materials. Another unique characteristic that 
distinguishes composites from other materials is that in the case 
of fiber composites, at the same time the structure or component 
is manufactured, the material itself is being synthesized. 



1.2 Stress Analysis for Composite Materials: 

A material as important as composite needs to be used with 
care. This need becomes all the more crucial, since the single 
most important reason for their usage is their lightness. Hence, 
any design which is oversafe, will increase the weight of the 
system disadvantageously . For better utilisation of the material, 
it is important that the stress field in the system is optimized. 
Hence the need for detailed stress analysis. 

There are several methods to find stresses in a composite 
laminate. Some of them are; photo-elasticity, strain-gage 
techniques, holography, moire’s method, FEM. 


1.3 Photoelasticity and Composites: 

Photoelasticity, a well established branch of experimental 
mechanics, has been successfully tried for the analysis of 
isotropic materials. 

The advantage of photoelasticity is that it is a whole-field 
technique. Consequently, one single experiment in photoelasticity 
gives a wealth of data. Hence the reliability of each experiment 
in photoelasticity is definitely much more than an experiment 
based on single-point techniques. 

However, isotropic transmission photoelasticity runs into 
trouble when composite materials are the field of study. It is so 
since composite materials are anisotropic in nature. 
Consequently, the laws of conventional transmission 
photoelasticity fail to predict the mechanical behaviour of 


composite materials. 
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In of "these facts, it becomes "very essential to modify 

the laws of isotropic photoelasticity. In this field, several 
investigators have made their significant contributions. 
Consequently, a newer version of photo-elasticity, "orthotropic 
photoelasticity" has emerged. 

Investigations in the field of orthotropic photoelasticity 
can be categorized into three broad areas. 

1. Development of suitable composite model material for 
photoelastic investigations. 

2. Development of suitable stress-optic laws and strain- 
optic laws which are accurate as well as easy to handle. 

3. Using the developments mentioned in 1 and 2 to analyze 
physically realizable problems, 

1.4 Fracture in Composites: 

The importance of fracture mechanics, especially linear- 
elastic-fracture-mechanics is well documented. This has been 
dealt in great detail by Liebowitz [9], Broek [103 » Kobayashi 
[11] and Rolfe [12], Although most of the work done in this area 
has been done for isotropic materials, the significance of 
anisotropic materials like composites in the modern industry, has 
impelled several investigators to establish the theory of 
anisotropic fracture mechanics, and furthermore, to investigate 
their fracture mechanisms. 

In these studies the model material is assumed to be 
homogeneous and orthotropic if either the structure is free from 
flaws which may be the cause of an eventual failure initiation, 
or the structure may have a flaw whose size is large in 
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comparison with the local micro- structural parameters such as 
fiber diameter and the distance between two adjacent fibers. 
Otherwise, in failure initiation studies the material has to be 
treated as a non-homogeneous continuum containing local flaws 
with certain geometries. 

A distinguishing feature of composites, from the view of 
structural failure, is that these materials are generally not 
isotropic with respect to their fracture resistance. Furthermore, 
in a majority of cases, the planes of orthotropy are generally 
also the planes of weak fracture resistance. Thus in orthotropic 
materials regardless of the overall geometry and the loading 
conditions, the fracture propagation would be either along a 
plane of orthotropy or would have zig-zag path. This is in sharp 
contrast to isotropic materials, where the plane of fracture 
quite often coincides with the plane(s) of maximum principal 
stresses. This is one very important reason for carrying out an 
in-depth study of the fracture-phenomena of composites. 

A very important parameter in fracture studies is the 
stress- intensity factor. Knowledge of this factor aids in 
predicting the stress-field near the crack-tip in the material 
under observation. However, usage of this parameter is justified 
as long as the zone of plasticity in the vicinity of the crack 
tip is small enough. For ductile metals, the zone of plasticity 
is very large. However, plastic yielding in the case of composite 
materials is almost negligible. Thus, the relevance of stress- 
intensity factor for the understanding of fracture in composites 
is all the more important. 



1.5 Present Work and Its Scope: 

The significance of stress-intensity factor for composites 
has been discussed in the preceding section. So far the method of 
photo-elasticity, despite its advantages as discussed above, has 
not been employed with much of success for the determination of 
this important parameter. 

One very important reason for this is that the availability 
of a suitable orthotropic material with a high degree of 
transparency is very difficult. As a part of the present 
investigation, a suitable photoelastic model material for 
orthotropic materials has been developed. The method though 
tedious, gives very good results. The details of this development 
are given in the second chapter of this thesis. 

In view of the tedious process of the fabrication of the 
model material, optimum usage of this material, especially during 
its characterisation, is very important. The relevant details of 
the characterization of the model material have been discussed in 
the third chapter. 

In the fourth chapter, stress intensity factor for 
unidirectional laminates has been evaluated, using orthotropic 
photoelasticity. The method followed for this purpose has been 
found to be quite reliable. However, the process is time 
consuming. The fifth and the final chapter suggests various ways 
of automating this process and making it even more rigorous. 
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2 DEVELOPMENT OF TRANSPARENT COMPOSITES : 

2 . 1 Int-roducbion : 

The first step in applying transmission photoelasticity 
method to the stress analysis of orthotropic material involves 
the development of a suitable model material. The properties 
desirable of a model material for photo-orthotropic-elasticity 
applications are : 

1. The model material should be transparent and homogenous. 

2. Volume fraction of fibers should be high enough to 
ensure sufficient orthotropy in the model. However, fiber 
volume fraction should not be too high so as to reduce 
its sensitivity. 

3. The model material should be designed to exhibit 
adequate amount of birefringence. Use of a very thin 
model will result in inadequate response of the stressed 
model to polarized light. On the other hand, use of a 
very thick model causes appreciable variation of the 
stress-field in the thickness direction. This is 
significant during oblique incidence. 

The first two requirements were met by controlling the 
fabrication process and varying the constituents of the model 
while the third requirement was met by controlling its thickness. 

2.2 Review of The Past Development Transparent Composites: 

First attempt to prepare transparent orthotropic material 
was made by Horridge [9]. He matched the refractive indices of 
glass and polyester resin approximately to achieve the desired 
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results. Vinyl di-chlorosilane was used to improve the glass- 
resin bonding. Pih and Knight [2] used E-glass fibers for this 
purpose. They u.sed wet filament winding process. The entire 
process was carried out in vacuum to ensure removal of air 
bubbles. Sampson [3] also made attempts in similar directions. 
However, he used dry filament winding process. 

Dally and Prabhakaran [43 prepared large sheets of 
transparent composites by using E-glass fibers and Paraplex P- 
444H resin. Thirty percent styrene was added to resin to ensure a 
close match between the refractive indices of glass and resin. 
Dally et. al. characterized these sheets as sufficiently 
transparent. FMC Corporation [5] has developed transparent 
composites by using two different special formulations. Usage of 
silane coupling agents is advised in these formulations to ensure 
that laminates "retain there clarity upon aging [5]“. Pipes and 
Rose [6] applied tension in the fibers to ensure proper alignment 
of fibers. However, the model material developed by them was not 
homogenous in nature, since the density of the fibers was larger 
near the mid-plain than in the surrounding areas. Agarwal and 
Chaturvedi [7] developed highly transparent composite materials. 
Their approach involved the matching of refractive indices bf the 
glass fibers and matrix material as closely as possible. This was 
achieved by using two additives; di -vinyl benzene and di -butyl 
pthalate to alter the refractive index of polyester resin to the 
desired level. Fairly good results were achieved in this way . 

The model materials developed by all of the above 
investigators bad been successfully used in their investigations. 
However, it is desirable that isoclinics and isochromatics 
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obtained from such models be still sharper for more accurate 
analysis. The problem becomes even more acute for dynamic 
fracture analysis. It is for this reason that an improvement in 
the optical properties of model material was undertaken. Such an 
objective can only be met by understanding the behaviour of light 
as it passes through a non-homogenous media, such as a composite 
laminate. 

2.3 Optical Considerations for a Model Material*. 

In section 2.1, various properties required for a good 
orthotropic material have been discussed. Of these, the first two 
conditions can be satisfied by controlling various parameters 
related to material processing. What follows are details of how 
various requirements for arriving suitable model material can be 
met. 

2.3.1 Transparency : 

The transparency of a model material depends on: 

1. Transparency of matrix material: The matrix material 
should be transparent. Fortunately, most of the epoxy and 
polyester resins, which are transparent and are been 
extensively used for photoelastic analysis of isotropic 
materials, can also be used as matrix materials. 

2. Matching of Refractive Indices of Matrix and Fibers: The 
refractive indices of matrix material and glass fibers 
affect the transparency of the composite signif 
Severe optical inhomogeneity (caused due to mismatch 
between refractive index of glass fiber and matrix 
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material) can render an otherwise good model material 
useless in two ways. Firstly, the loss , transparency 
occurs due to scattering of light at the fiber-matrix 
interface. Consequently, the model material becomes 
opaque or turgid. 

Another optical phenomena, which is more detrimental 
than the loss of transparency in model material for 
photoelastic application, is due to differential 
refraction of light passing through an optically 
inhomogeneous system. The effect of Snell’s Laws of 
Refraction may be illustrated by considering the path of 
light as follows: 

Consider that a parallel beam of light is incident on 
a body which has two parallel faces. The body may be: 

1. Optically homogeneous (Fig. 2.1-a). 

2. Having two different materials which have different 
refractive indices and an interface plane parallel 
to external surface on which light is incident (Fig. 
2.1-b). 

3. Having solid medixun with a solid inclusion which has 
a curved external surface (Fig. 2.1-c). 

In the first case, all rays of incident and emergent 
beam will be collinear and mutually parallel. In the 
second case, the angle of deviation for each ray in a 
particular medium will be the same. However, this angle 
of deviation will vary from medium to medium. Further, in 
this case, the emergent rays will be mutually parallel, 
subject to the restriction that incident rays are 
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parallel. In the third case, the situation is entirely 
different. Here, different rays will deviate by different 
angles. Hence, the emergent rays will not be parallel. 

In glass fiber reinforced materials, glass phase is 
uniformly dispersed in the matrix material in the form of 
thin circular cylinders. Hence, the cross-section of a 
composite plate containing a single glass fiber can be 
represented as depicted in Fig. 2.1-c. However, in actual 
systems, there are a number of such cylinders. 
Consequently, a beam of light emerging through such a 
composite plate will be of a highly dispersed and 
diffused nature. Thus, even a slight mismatch in the 
refractive indices of glass fiber and matrix reduces the 
suitability of composite plate for photoelastic purposes 
appreciably. Therefore, the first challenge to produce a 
superior photoelastic orthotropic model comprises of 
matching the refractive indices of glass fiber and matrix 
as closely as possible. 

3. Presence of Entrapped Air: Any air bubbles if present 
in the matrix phase, lead to same problems as does the 
mismatch between the refractive indices of glass fiber 
and matrix. Hence, it is very important to remove any air 
which gets entrapped in the laminate. One method to 
remove these bubbles is by curing the resin in vacuum. 

The affect of the presence of glass fibers in a composite 
model material on the path of light was studied in great detail 
by Agarwal and Chaturvedi [7]. They found out that: 
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1. A parallel beam of light appears smeared in the 
direction perpendicular to the fibers. 

2. Sharp boundaries of an object, when viewed through the 
model, appear smeared again in the direction 
perpendicular to the fibers. As a consequence, very 
closely spaced straight and parallel lines do not 
appear distinct. 

These two effects, the investigators concluded, cause a 
reduction in the sharpness of isoclinic and isochromatic fringes. 
Hence, an orthotropic photoelastic model is superior when it is; 

1 . Highly transparent . 

2. Causes minimum distortion of the view of an object when 
seen through it. 

2.3.2 Homogeneity and Orientation of Fibers: 

The properties of model material should be same at all the 
points in the body. Another requirement for a satisfactory model 
material is that the orientation of fibers should be as exact as 
desired. Both these conditions will also ensure the 
reproducibility of the results. 

2.4 Fabrication of the Model: 

As discussed earlier, the two most significant 
considerations for obtaining optically superior composite model 
materials constitute matching of the refractive indices of fiber 
and matrix and removing any air which gets entrapped. 

To achieve these ends, various combinations of matrix 
material and glass fibers were tried. 



13 


2.4.1 Selection of Glass Fiber: 

Among glass fibers; both, glass roving and unidirectional 
mat were tried. There specifications are: 

1. Unidirectional Glass Fiber Mat: 

Manufacturer : Fothergill Engineered Fabrics Ltd., U.S.A 
Style : 1200 Tex, 4 Threads per cm. 

Area Density : 500 gm,/mm 

Type : Unidirectional plain 

Width : 50 cm. 

2. Glass Roving: 

Manufacturer : Ceat Ltd. , Glass Fiber Division, Hyderabad 

Style : 300 Tex 

Brand : DR 600 FE 

Size : Each spool of 4.0 Kgs. 

Not much success was achieved with glass rovings. The 

principal reason for this was poor penetration of resin 
compounded with poor wetting of fibers. Problems encountered with 
the usage of unidirectional glass fiber mat were of a different 
nature. In this case, there was little problem as far as 
penetration and wetting was concerned. However, using glass fiber 
mat as it is, led to opacity in some zones of the laminate. This 
was due to the presence of thick pre-impregnated bundles of glass 
fibers in warp and weft directions at these locations. These 
fiber-bundles have been introduced in the mat by the 
manufacturers to maintain the orientation of glass fibers. 

In light of these observations, it was realized that glass 
fiber mat is preferable over glass roving. It was further 
realized that the model material developed from unidirectional 
mat will be perfect only if the “thick bundle of prepreg fibers" 
could in some way be removed from the mat. 
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2.4.2 Selection of Matrix System: 

Four grades of matrix materials were tried. They were: 

1. Epoxy resin {LY-5052) and hardener (HY- 505 2) manufactured 
by Hindustan Ciba Geigy Ltd. , Bombay . 

2. GP (general purpose) polyester resin manufactured by 
Parikh Chemicals Ltd. , Kanpur. 

3. Clear Casting Grade Polyester Resin manufactured by 
Parikh Chemicals Ltd. , Kanpur. 

4. Polyester resin having mono-maleic acid (MMA) monomer. 

Each of these resins was used with glass roving as well as 

glass fiber unidirectional mat. In addition, suitable silane 
coupling agents were also used in these matrix - systems . The 
influence of s.ilane coupling agents was not dramatic. However, 
wetting improved to some extent in cases where glass roving was 
used. 

Wetting and penetration were good enough in case of epoxy 
resin. However, problems encountered with this resin were: 

1. The colour of matrix material after curing was dark 
yellow; which is not preferable. 

2. The time required for the resin to solidify was not 
large. A longer curing time is helpful in removal of 
entrapped air bubbles, when casting occurs in vacuum. 

The advantages of clear casting grade polyester resin were 
its low viscosity and colorless complexion. However, the 
difference in the refractive indices of fiber and matrix was 
found to be so large that it was not considered worthwhile to 
work with this resin system. Similar problems were also 
encountered when polyester resin having MMA polymer was used. 
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General purpose polyester resin was found to be the most 
suitable resin. Its refractive index was very close to that of 
glass fiber. The small difference between refractive indices of 
fiber and matrix was further minimized by adding 10% styrene by 
weight . 

2.4.3 Fabrication Procedure for the Model Material: 

Initially the ends of the fibers were bound using epoxy 
resin. For this purpose, epoxy resin (LY556) and hardener 
(HY951), manufactured by Hindustan Ciba Geigy Ltd. were used in 
the ratio of 10:1 respectively. Such a pre-treatment process was 
"essential", since in the absence of this treatment, the 
alignment of the fibers gets disturbed while treating the mat for 
removal of pre-preg fiber bundles. The relevant details of pre- 
treatment procedure are shown in Fig, 2.2. 

Next, the matrix material was prepared. For this purpose, 6P 
polyester resin was used. This resin has a low viscosity. Methyl 
ethyl ketone peroxide (MEKP) was added as 1.0 percent of the 
weight of resin to act as catalyst. 1.0 percent cobalt octate (by 
weight) was used as accelerator. The resin was heated in a 
temperature controlled oven to ensure the removal of traces of 
water droplets, present if any. It was also ensured that 
accelerator was added to resin after the addition of catalyst. 
These two chemicals were mixed uniformly in the resin. 

Pre-treated unidirectional glass mat was cut in pieces of 
desired size. Next, each piece of mat was soaked in resin. Care 
was taken to minimize the amount of air entrapped at this stage. 
Hence, resin was spread on the mat with a rubber lined roller, 
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Fig. 2.2 Zone in which epoxy resin is to be applied 
for Pre-Treatment process. 
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instead of a brush. Once the mat was fully wet, bundles of pre- 
impregnated fibers were removed using a pair of tweezers. Removal 
of pre-preg fiber bundles is possible only when the fibers are in 
a wet state, since the adhesive which binds pre-preg fibers to 
the mat, gets dissolved in polyester resin or styrene or acetone. 
A parallel attempt was made to remove these pre-preg fibers and 
then applying acetone or styrene to the mat, followed by drying 
the mat in oven, and then using the mat for casting purposes. 
However, not much success could be attained in this direction 
since the application of styrene or acetone led to the removal of 
'size’ from the surface of glass fibers and hence caused poor 
wetting. 

In this way, several layers of glass fibers were stacked up, 
and covered with mylar sheet on both sides. The entire set up was 
placed between two perspex plates separated by 3 mm. thick 
aluminium spacers, as shown in Fig. 2.3. 

Next, the fibers were subjected to tension to ensure proper 
alignment of the fibers. This was done by clamping the ends of 
the fibers by a steel plate and adjusting the position of the 
tension screw (Fig. 2.3). This was followed by clamping the 
perspex plates by nuts and bolts. An interval of 5-10 minutes was 
given to facilitate the flow of extra resin out of the mold, due 
to clamping pressure. Later, the tension in the fibers was 
removed and the entire setup was subjected to vacuum pressure for 
approximately five hours to ensure the removal of air bubbles. 
This was followed by curing the resin at 60 for 12 hours 
followed by postcuring the laminate at 80 for six hours. 




Fig. 2.3 Diagram of cas-bing p 
to subject fibers to tension. 
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2.5 Optical Properties of Model Material: 
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Or-tical properties of uni directional ooiaposite model 
jnateriiil was e;:>.ii»iined to assess the suitability for rhotoelastic 
c.i::i3yoia .Mid allied applications. The ruodel material has been 
t'-.-tc'3 f<<r transi>arerioj> , resolution of closely spaced lines and 
pho^ Oil Last ic effect. The relevant details follow: 


2.5.1 Tranaparenoy : 

Tl^e transparency of the W'^del material developed was 
a;-', -..cs.-ried in toms of transmission ratio. Transmission ratio is 
defined as tho ratio of incident light’s inten.jity to that of 
•-mereent light. 

Tho circular disc model was placed in the path He-He baser 
l The omorgisnt beam was focussed, using a convex lens , on 

th^ sensing element of the photometer. Inten.sities of the 
em':.rt,cnt ouU incident bean* vrere measured using a digital 
photometer. The entire setup is shown in Fig. 2.4. The value of 
oraiismission ratio for the model material (3.2 rom thick) 
developed was found to be 0.69 and was greater than that for tlie 
modr-1 material developed by Chaturvedl [8] which bad a 
t.ransjjil sj,ion ratio of 0.53. 

2.5.2 Eer.olution of Lines Through the Model Material: 

J t ha.s boen discussed in the preceding passages that the 
sharp boundaries of an object appear blurred when viewed through 
a composite model whose constituents do not have well matched 
refractive indices. This effect was exaiftined by viewing closely 
spaced straight and parallel lines through the models. Two st»ts 
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of equispaced lines were drawn inside a circle with lines in one 
being perpendicular to those in the other set (Fig. 2.5-a). These 
lines were viewed through through circular discs of 
unidirectional model material (Fig. 2.5-b). Direction of the 
fibers has been indicated by arrows in the figure. It is clear 
from the figures that the lines are distinguishable from one 
another. Placing an optically inferior model over such a set of 
lines revealed that lines which were parallel to the fibers 
appeared blurred, 

2.5.3 Photoelastic Effects: 

The model materials have been examined for photoelastic 
effect by placing their circular disc models in a circular 
polariscope employing a diffused sodium light source. The discs 
were loaded in diametral compression with fibers oriented at 90° 
to the loading axis. The resulting fringe patterns for 
unidirectional models have been shown in Figs. 2.6-a and 2.6-b, 
the latter pertaining to an inferior model material. 

It is amply clear from the figure that the fringes for 
superior model material are much better defined than for other 
material (Fig. 2.6~b). Some nonuniformity of fringes parallel to 
fibers occurs due to non uniform distribution of fibers. For this 
reason the photoelastic sensitivity of the model material varies 


from point to point. 
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Fig. 2.5: View of parallel lines through a; 

(a) Superior model material. 

(b) Model with larger mismatch of refractive indices. 
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Fig. 2.6: Isochromatic fringe patterns for; 

(a) Superior model material. 

(b) Model with large mismatch of refractive indices. 
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3 CC»1P0SITE CHARACTERIZATION : 

3 . 1 Introduction : 

Any material, when newly developed needs to be fully and 
accurately characterized before it can be put to other uses. This 
necessity becomes all the more important for composite materials, 
since their properties depend on a large number of parameters. In 
the present investigation, the newly developed material was 
tested for its mechanical and photoelastic properties. The 
former included elastic constants and Poisson’s ratio and elastic 
constants, while the latter included three stress-optic 
coefficients . 

3.2 Experimental Setup : 

While fabricating the model material it was felt that the 
process of fabrication is very tedious and time consuming. For 
this reason, it was natural to design and develop a suitable 
experimental setup to ensure: 

1. That the material be characterized accurately. 

2. The amount of model material consumed is minimum. 

The available infrastructure in the laboratory included a 
MTS machine (suitable for mechanical characterization) and a 
polariscope oxm loading frame (suitable for photoelastic 
characterization). One option available, was to use these two 
facilities separately for their respective purposes. However, 
this was later found to be unfeasible due to the following 


reasons: 
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1 . The material developed seemed to a have very high value 
of stress optic coefficient in the longitudinal 
direction (FqL) . This in turn meant that very large 
loads were needed to produce appreciable photoelastic 
effect in the material. However, the conventional 
loading frame available was not found to be strong 
enough to bear such high loads. 

2. A procedure for characterising the model material 
separately for its photoelastic and mechanical 
properties required larger amount of model material an 
time. 

In light of the above constraints, it was decided to develop 
a new setup which could suffice the following conditions: 

1. It was rigid enough to take high loads. 

2. It incorporated certain features which could facilitate 
photoelastic as well as mechanical characterization of 
the model material simultaneously. 

The new setup comprised of the MTS machine having a quarter 
wave plate and a polarizer on its one side and an analyzer along 
with another quarter wave plate on its other side. A light box 
having white and sodium light sources was kept behind the 
polarizing element. While arranging the polariscope’ s elements on 
either side of the MTS machine it was ensured that surface of the 
tension test specimen when loaded on the MTS remained parallel to 
the optical elements and also that the light coming from the 
light box was incident normally on the specimen. Such an 
alignment between the specimen, elements of polariscope and light 
box was maintained by using a try-square. This was very important 



since laws of photoelasticity are valid as long as the angle of 
deviation of a light ray passing through the specimen is sero. 

3,2.1 Material Testing System: 

The Material Testing System (model MTS 810) is a servo- 
hydraulic testing system designed with a view to accommodate a 
wide variety of standard material tests. The testing machine has 
the load range capacity of -100 kN to +100 kN which can he 
applied in the form of a sine, ramp, haversine, triangular and 
other programmable waveforms. It uses a microconsole (458.20) to 
provide closed loop control of the servo-hydraulic system. 
Specimen is gripped in the load unit. Hydraulic power supply is 
linked to the microconsole. To control the tests, microconsole 
has a Displacement AC controller. Load DC controller and Strain 
DC controller. Testing can be done in any of these control 
modes. To generate a waveform microconsole, has a microprofiler 
which can be operated in three modes viz. programmed, direct and 
remote. This machine can perform Tensile, Fatigue and some other 
type of tests. 

The entire setup has been schematically illustrated in Fig. 
3.1. A photograph of the setup is shown in Fig. 3.2. 

With the aid of such a setup, it was possible to record the 
following information during the tension test: 

The external load on the specimen. This was 
accomplished by using a load-cell having a maximum load 
rating of 10 kN. 


1 . 
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Fig. 3.1 Schematic layout of experimental setup 
developed for material characterization. 
















Fig. 3.2: Photograph of experimental setup prepared 
for material characterization. 
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2. The strain at any desired point on the specimen. This 
could be found out by using an extensometer meant for 
measuring longitudinal strains. 

3. The photoelastic response of the material in terms of 
the fringe order at the point of interest. For 
recording the fringe order, monochromatic sodium light 
was used. White light was not preferred since the 
ordering of fringes for such a light is very difficult. 
This is so, since the colour of fringes changes 
appreciably with respect to change in fringe order 
Hence, recognition of these colours depends heavily on 
human judgment, which can introduce very large errors. 

From this information, extensional elastic moduli and 
stress-optic constants for any particular direction could be 
evaluated simultaneously. 

By testing samples with fiber orientation at 0*^,90'^ and 46°, 
the values of E^, Et, E 45 , F^,, F£,,F 45 , were found out. However 
for the evaluation of shear modulus (Glt) and Poisson’s ratio 
(VLT) > lateral strains were needed. These strains were measured 
by using lateral extensometer. The setup of lateral extensometer 
is shown in Fig. 3.3. 

3.3 Hechanical Characterization : 

The volume fraction of fiber the for model material has been 
calculated as follows: 

Thickness of the specimen = 3.2 ram. 

Density of glass fiber = 2.54 gm./cc 




32 


Area density of the mat = 405.0 gm. per sq. m 
No. of layers in each laminate = 7 

For 1 m® of laminate, volume of glass= 405-7/2.54 = 1116.1 cc 
Volume of laminate = 3200 cc. 

Volume fraction of fiber = 1116.1/3200 = 34.88 % 

Other properties which were required for the model material 
are : 

1. Extensional modulus in the longitudinal direction (El). 

2. Extensional modulus in the transverse direction (Ef). 

3. Poisson’s ratio (vlt)« 

4. Shear modulus (Glt). 

For evaluation of these properties the following relations 
can be used: 


^X - “^x/^x 

when 

o 

it 

Txy = 0 

(3.1) 

"'^xy “ “^y/^x 

when 

O-y = 0, 

"^xy ~ ® 

(3.2) 

Glt “ ^ ^ ^"^^xy ^ ^ 

when 

CTy = 0 , 

'■'xy = 0. cfx ~ *^45 

(3.3) 


Using Eq. (3.1), elastic moduli (El, Et) can be easily 
evaluated. Eqs. (3.2) and (3.3) can be employed to find out the 
Poisson’s ratio (vlt) shear modulus (Glt) respectively. 

For the calibration of the model material, static uniaxial 
tension tests were carried out on flat specimens with rectangular 
cross- sections . V?hile conducting the tests, the strains in the 
longitudinal and lateral directions were recorded along with the 
applied loads. The details of the experimental setup have been 
given in section 3.2. 

By conducting uniaxial tension tests for 0*^ and 90*^ samples, 
El, Et and vlt were evaluated. For finding Glt. off-axis 45^ 
coupon test was conducted, and using Eq, (3.3) Glt was computed. 
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While conducting the test in the longitudinal direction, it 
was ensured that the direction of application of load coincided 
with the material direction of the orientation of the fibers. 
This is very important since the gradient of elastic modulus in 6 
direction (Eg) with respect to orientation of external load (i.e. 
dEg/dS) is very high near 0*^, 

3 - 3 Photoelastic Characterization : 

3.3.1 Assessment of Existing Theories of Photo-Orthotropic 
Elasticity: 

A polarized study of glass fiber laminates by Horridge [9] 
seems to be the first attempt in the direction of photoelastic 
application to reinforced plastics. Pih and Knight [2], proposed 
a stress optic law for transparent composites. They also observed 
that isoclinics are not representatives of the principal stress ; 
directions. 

1 

Fere = o • t/Ng I 

However their stress-optic law, based on the concept of 
average optical behaviour and involving a single photoelastic 
constant for orthotropic materials, was unrealistic. 

Sampson [3], using the analogy of the Mohr’s circle for 
stress suggested an equivalent circle for birefringence. He 
postulated that orthotropic materials can be characterized with 

I 

! 

three principal stress. values Fq-l, Fcyf, Fq-lx^ ^ind that the stress i 
system ert,' ”^LT produces three components of birefringence 

which are defined as: 
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Nl = ctl/^oL (3.4-a) 

NT - oy/FaT (3. 4 -^ 1 ) 

NlT = S^tlt/FcjLT (3.4-c) 

where the subscripts L and T refer to principal material 
axes. The resultant birefringence is: 

N = [{ctl/Fol ~ o'T/^’aT)' (3.6) 


Dally and Prabhakaran [4] refined the stress proportioning 
method used by Pih and Knight [23, to derive another stress-optic 
law. Their formulation was governed by the following equation: 

^ - Njjiatrix ^fibers (3.6) 

From Eq. (3.6), the authors established equivalence between 

the values of Fe as predicted by their theory and that by Eq. 
(3.5). This is given in Eq. (3.7). 

(^ePrabhakaran)®* = (5'eSampson)® (3.7) 

where 0 is a difference term. 

It was found by these investigators that the maximum value 
was 1.06 when 6 assumed a value of 28 deg. Thus if the 
difference term in Eq. (3.7) is neglected, since it is very 
small, then then the formulation of Sampson’s and Prabhakaran’ s 
are identical. However, these investigators did not discuss the 
physical significance of the isoclinics. 

Bert [143 assigned tensorial nature to birefringence on the 
basis of crystal theory of Bhagavantham [153- By relating a 
second order birefringence tensor to second order stress tensor 
through a fourth order photoelastic-property tensor, Bert 
formulated a stress-optic law. However, his relations were not 
simple enough to be used in experiments. Besides this, Bert did 
not corroborate his predictions by experimental evidence. 



Pipes and Rose [6] proposed that strain-optic laws could 
only be employed for composites with low volume fraction of 
fibers and verified their contention experimentally. Further, 
they developed a strain-optic law by assuming that birefringent 
composite materials obey the classical strain-optic relations for 
isotropic materials. For unidirectional fiber composites , the 
strain-optic law proposed by them was: 

- €2 = n*f€/t = N-Fc (3.8) 

In Eq. (3-8), subscripts 1 and 2 indicate the direction of 
principal strains and Fg is the strain-optic coefficient. A novel 
feature of their law is that it requires only one material 
constant. However this relation can not be used for composites 
having large fiber content. 

Agarwal and Chaturvedi [7] have assessed various theories of 
photo-orthotropic-elasticity on the basis of an extensive work 
conducted on 'superior* model materials developed by themselves. 
They have shown that the isoclinics are closer to the contours of 
principal strain differences rather than those for principal 
stress differences. This very important observation prompted 
Agarwal and Chaturvedi [16 3, to derive a general strain-optic law. 
This derivation was mathematically sound. Their exact strain- 
optic law based on Bhagvantham’ s [15] theory of crystalline 
photoelasticity is: 

N = [{^L/^€L “ ^T/^€T5' { 2‘5TjT/^€LT3- 3^ 

Tan2v6 = [yLT/^€LT3 C^L/^eL " ^t/^ 6T3 
where ; 

N = Fringe order 

i> - Isoclinic angle 


(3.9) 

(3.10) 
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F = Strain-optic coefficient 

y - Shear direction 

€ = Extensional strain 

L,T K Refer to longitudinal and transverse directions 

Further, these investigators realised that the degree of 
orthotropy with respect to strain-fringe value is considerably 
smaller than with x'espect to stress-fringe value. This finding 

prompted the investigators to conclude that an approximate 

strain-optic law was likely to introduce much smaller errors than 
an approximate stress-optic law. On the basis of this intuition, 
they proposed an approximate law to estimate the direction of 
principal strains and also to predict the values of principal 
strain difference at a point. It was found that the error 
introduced by such an approximate law was small enough and within 
the experimental limits. 

r 

N = C{€l - € 1 -} + {2-yL^^ -i- F^ (3.11) 

Tan29i = y^T (3.12) 

Agarwal [17] proposed another approximate strain-optic law 
which was much more accurate than as given in Eqs. (3.11) and 


(3.12). This law was 


obtaine 


approximating 


ratio (FgL/Fex) unity on the basis of experimental evidence. 
The resulting modified law becomes: 

N •- (€l - € 7 ) • [cos® 26 + (F€L/F€LT)®“si 3 a® 2 e]’^ -J- F^ (3.13) 
Tan20 = (Fl/Flt) * Tan26 (3.14) 
With so many theories available, each varying in its 


accuracy and approachability, it required sound judgment to 


choose the most appropriate one of these for further 


investigations. Two considerations which governed the selection 
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of such a theory were: 

1 . Accuracy of the theory . 

2. Ease in usage of the selected photo-orthotropic law. 

In light of these two constraixxts, it was realized that the 
determination of three stress-optic constants (F^l, Fq-lt) as 

proposed by Sampson for the model material would suffice our 


purpose. These three constants are defined as under: 

FoL = (3.15) 
F 0 T - "(3.16) 
Fae = oe-t/Ne (s.i?) 
FoLT ^ - 0.25-(l/fT - l/fL)^3“^ (3.18) 


Thus, using Eqs. (3.15) and (3,16), Fql FoT were 
evaluated, while Eqs. (3.17) and (3.18) were employed to find out 
FqLT evaluated. The exact fringe order at any point was found 
out by Tardy’s method of compensation. The details of the 
experimental setup are given in the section that follows. 

3.5 Analysis of the Collected Data: 

3.5.1 Evaluation of Extensional Elastic Constants and Stress 
Optic Coefficients: 

Currently, most of the investigators characterize a material 
with point-matching techniques. In such techniques it is assximed 
that materials behave in exactly the same way as predicted by 
mathematical models and thus the errors involved in the 
experimentation are not ^arge. Thus very few data points are 
collected for analysis. Very little significance is attached to 
the amount of errors which may creep up in such an analysis. 
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These errors may be due t-o residual stresses, residual strains or 
various personal factors. This problem becomes all the more 
significant while calibrating the material for its photoelastic 
properties. In the present investigation, the method of least- 
square analysis was employed to filter out the error introduced 
due to various reasons during such characterization. 

3 . 5 . 1 . 1 Linear Least Square Method : 

Consider the following relations: 

Oi = Ei-€i (3.19) 

Oi = (l/Fcyi)‘(Ni/t) (3.20) 

where; 

= Externally applied stress along direction i 
F^yi = Stress-optic coefficient for direction i 
E^ = Extensional elastic modulus for direction i 
= Fringe order in the specimen for direction i 
= Strain in i direction due to o'i 
Eqs. (3.19) and (3.20) can be generalized as: 

S = M-R (3.21) 

where; 

S = External stress 

M = Material property 

R = Response of the material (i.e. strain or fringe order) 

If an additional constant term is added to the R.H.S of Sq. 

(3.21) to account for errors, then the said equation can be 
rewritten as: 

g=y-a«x+b (3.22) 

For a perfect fit, g is zero. However, in other cases 
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Eg.. (3.22) can be expressed in bbe mabrix form as: 

{g} = [c3-{s} (3.23) 

where ; 

{g} = Column vector for errors 

[c] = Rectangular matrix having individual data points 
{z} = Column vector for a and b as given in Eq. (3.11) 

Using Eq. (3.23), the values of a and b can be evaluated as 
follows: 

[c]'^-{g} = Cd]{s} 
where ; 

Cd] = Square matrix [cl'^-Eo] 

Thus ; 

{2> = [d] ^•[g3'^*{2> 

Using such an approach, extensional elastic moduli and 
photoelastic constants were evaluated. The values obtained were 
quite consistent with those mentioned in the available 
literature. The plots between stress aiid strain and between 
stress and fringe number for 0*^, 90° and 45° are shown in Figs. 
3.4, 3.5 and 3.6 respectively. 

3.5.2 Determination of Poisson's Ratio and Shear Modulus: 

Using a lateral extensometer a plot between lateral 
displacement and external load, was obtained for 0° and 45° 
specimens. From this plot another plot between lateral strain and 
external load, as shown in Fig. 3.7, was obtained. Next, a plot 
between longitudinal strain and load was obtained using a 
longitudinal extensometer. From these two plots and Eqs. (3.2) 
and (3.3), the values of vlt and Gli’ were computed. 
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3 . 6 Results : 

The values of various mechanical and photoelastic 
for the new model material were found to be: 

Vf = 34.88% 

El = 25.292 GPa 
Et = 4.6534 GPa 
E45 =: 5.4741 GPa 
vlt - 0.3877 
v-TL = 0.0712 
V45 = 0.3612 GPa 
Glt = 2.1103 GPa 
Fl = 143.67 N/mm 
F-f = 66.591 N/mm 
Ftl ~ 54.067 N/mm 


constants 




o 













Fibre 



o 


strain 
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4 KVALDATIOS OF MODE I STRESS INTENSITY FACTOR: 

4 . 1 Int.roduct.i on : 

The subject of fracture-mechanics is concerned with "the 
study of fracture of materials by crack propagation. The 
fracture-mechanics approach has been found to be very useful in 
predicting the strength of isotropic and homogeneous materials in 
the presence of flaws or cracks and also in developing techniques 
for quality control and in-service inspection. Fracture mechanics 
when applied to composite materials is much more complicated 
because self-similar crack growth, which occurs in homogeneous 
materials, usually does not occur in composite materials. Another 
factor which complicates the fracture-mechanics of composites is 
the presence of local heterogeneity in the material. 

It has been has shown that the strain-energy release rate 
for a material during crack growth can be correlated to the the 
stress distribution in the neighborhood of the crack tip [11]. 
Further, for isotropic materials, the stress field near a crack 
tip can be completely defined in terms of a parameter termed, 
stress intensity factor. This parameter is a function of the 
geometry of the specimen and external loading conditions. Thus, 
experimental determination of this parameter is very useful in 
understanding the phenomena of fracture. 

Unlike in isotropic materials, for anisotropic materials, 
the elastic constants vary along with direction. Hence, for 
composites, the fields for stress, strain and displacement in the 
vicinity of the crack tip are much more complex. 
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4.2 Orthotropic Linear Elastic Fracture Hechanics: 

The stress analysis ahead of a crack tip in composite 
materials is a very difficult task as compared to that for 
isotropic materials. The reasons for this have already been 
discussed in section 4.1. 

In view of these constraints, arriving at a closed form 
solution for stress field in the vicinity of a crack tip for 
composites materials is possible only, if it is assumed: 

1. The material is homogeneous and orthotropic. This 

assumption may be justified on the ground that the 
material may have a flaw (the crack length) whose size is 
large as compared to that of the local microstructural 
parameters such as the fiber diameter and the distance 
between two adjacent fibers. 

2. The crack growth is self similar in nature. 

With these assumptions, Irwin [183 and Wu [193 have employed 
Westergaard’ s differential equation [203 to find the stress field 
close to the crack tip. The governing equation is: 

+ 3%l/axi^X2^ + 8%1/^X2^ = 0 

Here is the Airy’s stress function. 

In this derivation the roots of the auxiliary equation: 

- b2*s^ bj^ = 0 (4.1) 

are required. Here: 

bi = ai 1/322 
b2 = [2-ai2 + a663/a22 
ail = 1/Ell 

ai 2 = “V12/E11 = -V21/E22 
322 = I/E 22 
aee = 1/^12 


(4.2) 
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The roots of Eq. (4.1) are either complex or purely 
imaginary. It was found that for our model material having 
elastic properties as mentioned in section 3,6, the roots of Eq, 
(4.1) are purely imaginary. Thus these roots may be expressed as: 


Si = i(3i S2 = -iPj^ - 

si = i|32 S2 = -i|32 - 

With this information Wu [19] derived the crack-tip 
equations as: 


(4.3) 

stress 


Ox = {Ki«|3i02*t32*cos(Y2/2)/'fE2 ” l3i'cos(Yi/2)/4Ri]/4( 27tr) + 

•cos(Yi/2)/^Rl - 32^’sin(Y2/2)/4R23}/(3i “32) + 
Cfox {4.4-a) 


oy = {Ki-{;3i*cos(Yi/2)/4Ri - 32‘Cos(Y2/2)/4R2]/4(2xr) + 

Kii-[sin(Y2/2)/4R2 - sin(Yi/2)/4Ri3}/(3i -32) (4.4-b) 

T^xy = Kii*[32'cos(Y 2/2)/4R2 - 3r cos(Yi/2)/4Ri]}/(3i -32) + 
{Ki-3i32'Csin(Yi/2)/4Ri - sin(Y2/2)/4R2]/4(2xr) (4.4-c) 


where ; 

Rj =. 4[cos®8 + 32*sia®0] (4.5) 
Tan Yj = 3j‘Tane (4.6) 
with 0 = 1 , 2 . 

In Eq. (4.4-a) the a^x was added by Mootahed and 
Zachary [21]. This additional term takes care of the effect of 
far field stress applied parallel to the crack. Further, if the 
crack is parallel to one of the material axes, then the terms 
involving Kn in Eqs. (4.4-a), (4.4-b) and (4.4-c) are irrelevant 
since in such a case. Mode II is absent. 

Eqs. (4.4-a) to (4.4-c) suggest that for the determination 
of stress intensity factor (S.I.F) one roust know the material 
propei'ties (3i, 32) the state of stress near the crack tip. 
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The former can be evaluated from the material properties given in 
section 3.6 of this thesis, while the latter can be evaluated by- 
employing various theoretical or experimental techniques. In the 
present investigations, the stress field was found out by- 
employing the principles of orthotropic photoelasticity. 

4.3 Method for Evaluation of Ki: 

4.^.1 Synthesis of Stress Optic Law and Wu’s Formulation: 

Sampson’s stress-optic law is a convenient tool for finding 
the stress field in the vicinity of a crack tip. The problem of 
residual stresses, which had not been addressed by Sampson has 
been taken care of in the following equation: 

(N/t)® = [oL/FaL -ctt/^oT + Nr*cos28r/t3® + 

[2«XL<r/^oLT Nj.*sin28i*/t]® (4.7) 

where ; 

N = fringe order 

t = model thickness 

= residual fringe order 
By = residual isoclinic angle 
= stress-optic coefficient 

By incorporating Eqs. (4.5) and (4.6) inEq. (4.7) fringe 
order (N) can be expressed as: 

(N/t)® = g(Ki, Oox' ^1’ ^ari» ®r) (4.8) 

For a given model (3i, |32> are known constants. Thus, Eq. 

(4.8) can be rewritten as: 

(N/t)“ = g(Ki, Oox» 6r) 
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Thus the fringe order at any point depends on Kj , r, 8, 
Nj- and 9j-. Hence, by finding the fringe order at a large number 
of points near the crack tip and by employing Eq. (4.9), Ki, Oqx* 
Nj. and 0j- can be found out. The numerical technique used for this 
purpose is the method of non-linear least squares [22] . 


4.3.2 Non-Linear Least Squares Method^ 

Eq. (4.9) can be rewritten as: 

f(Ki, Oox* 6< N/t, Nr, 8^-) = 0 (4.10) 

In this equation r, 8, N/t are known while Kj , crox> ©r 

are to be determined. The method of non-linear least squares 
essentially consists of : 

1. Linearizing Eq. (4.10) by a Taylor’s series expansion. 

2. Applying the criteria of least square to the linearized 
equations and finally iterating to converge the desired 
results . 

Eq. (4.10) can be rewritten in terms of unknowns as: 
f35(Ki, Oox’ Nr, ® (4.11) 

where k = 1,2, (M > 4) 

Taking Taylor’s series expansion of (Eq. 4.11) yields: 


(fk)i+l = (fk)i + ( ^^k/ SKi)i-hK]- + ( afk/5crox)i*horox + 

Ofk/^Nr)i-ANr + Ofk/^er )i'A035 (4.12) 

where; i refers to the i^^ iteration step and Kj, Oox> 

Nr, ©r corrections to their previous estimates. For the best 

fit, (fkH+1-*^- Eq. (4.12) yields: 


-(fk)i 


my ^Ki)i-AKi + 

(afk/aNr)rA^r + 


‘cl?1vrR 

I ! T 

dec. Nn A 


fL Lf’^RARY 


(4.13) 



Eq. (4.13) can be written in matrix form as: 


{f> = [b]‘{Ac} ^414 

where ; 

{f} = Error vector represented by the RHS of Eq. (4.13) 

[b] = Rectangular matrix containing partial derivative 

terms of Eq. (4.13) 

{Ac} = Correction vector containing Kj, etc. 

The correction vector {Ac} can be evaluated from Eq. (4.14) 

as : 

{Ac} = [d3'^*[b3'^*{f} (4.15) 

where ; 

[d] = [b3’^-Cb3 

The algorithm for evaluating the four unknown parameters in 
Eq. (4.11) is : 

1. Assuming initial values for Kj, Cox^ ^r> ®r- 

2. Computing the elements of matrices {f} and [b} for each 
of the data point. 

3. Computing {Ac} from Eq. (4.15), 

4. Revising the estimates of iinknowns i.e.; 


(Kl)i+l = 

(Ki)i + (AKi)i 

(4.16-a) 

(<^oxH+l ~ 


(4.16-b) 

(Ny)i+1 = 

(Nr)i + (ANr)i 

(4.16-c) 

(6r)i+i = 

( 9 J-) + ( A 6 j.) 

(4.16-d) 

Repeating 

steps (b), (c), 

(d) till {Ac} becomes 


acceptably small. 
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4.3,3 Initial Guess for Kj and Oox- 

The algorithm discussed in section 4.4.2 requires a good 
initial guess for the values of Kj and Oqx* initial guess 
is wide off the mark, then convergence to the exact value is 
extremely slow. Thus, the following procedure was adopted to get 
reasonably good estimates of Kj and Oox- 

Putting 6=90 and r=constant; Eqs. (4.4-a) to (4.4-c) give: 

CTjj^ = A-Kj + Oqx (4.17-a) 
Oy = B-Ki (4.17-b) 
Txy = 0 (4.17-c) 
where A and B are constants when r=constant. 

Eqs. (4.17-a) to (4.17-c) are simultaneous linear equations. 
They can be solved for Ki and Oqx "with aid of Eq. (4.7) by 
considering two points in the stress field lying on the contour 
6=90°. 


4 . ^ Experimental Details: 

Two types of specimens were used in this investigation. The 
geometry of the specimen is shown in Fig. 4.1. One set of model 
had fibers aligned parallel to the crack-length, while the other 
set of specimens had fibers running in a direction perpendicular 
to the crack length. The crack was simulated by using a 0.15 mm 
milling saw- cutter. Each crack was 6 mm long while the breadth of 
the specimen was 25 mm. Thus the ratio of crack length (a) 
specimen width (w) was 0.24. This was necessary since for larger 
values of a/w ratio, free edge effects modify the stress field 

signif icantly. 




The optical and elastic properties of the model material 
have been mentioned in section 3.6 of this report. For such a 
model, the values of nondiraensional material properties (3i. 02 
when crack is perpendicular to the fibers, are: 

01 = 1.40282 

02 = 0.30567 

Single edge notched (SEN) specimens were subjected to 
tensile stresses. This was done using the experimental setup 
described in section 3.2 of this thesis. 

The isochromatic fringe patterns for each specimen were 
recorded by using a photographic camera. These patterns are shown 
in Figs. (4.2), (4.3), (4.4). These photographs were later used 
to collect data points. From each photograph, the details of 30- 
35 data points were collected. The co-ordinates of each of these 
points and the corresponding fringe orders were used as the input 
data field for a program written in FORTRAN used to evaluate the 
mode I SIF for each of the SEN specimens. The initial value of Kj 
and Oox was calculated in the same way as explained in section 
4.1.3. The algorithm of the program was same as that given in 
section 4.3.2. A listing of this program is given in Appendix 1. 

4. S' Results: 

Following the procedure detailed in section 4.5, the 
following values of S.I.F in mode I were evaluated: 




Fig. 4.3: Isochroiaatics for SEN specimen Si 

(a) OL ” 37.2 MPa 

(b) CL = 49-6 MPa 


Fig. 4.4: Isochromatics 

(a) 0L = 37.2 MPa 

(b) OL = 49.6 MPa 


for SEN .specimen S 2 
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Fig. 4.5: Isochromatics for SEN specimen S 3 

(a) OL = 37.2 MPa 

(b) OL = 49.6 MPa 
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For 0 '^ specimens with following material properties: 
Vf = 34.88% 

El = 25.292 GPa 
Ef = 4.6534 GPa 
E 45 = 5.4741 GPa 
’^LT " 0.3877 
v-tl = 0.0712 
V 45 = 0.3612 GPa 
Glt = 2.1103 GPa 
Fl = 143.67 N/mm 
Ft = 66.591 N/mro 
Ftl ~ 54.067 N/mm 
Initial guess for Kj = 6.2 MPa4m 
Initial guess for =49.4 MPa 
the values of S.I.F and o^x were : 


Stress (ol^ 

S.I.F (Kj) 

Far Field Stress (o^x) 

(MPa) 

(MPa4m) 

(MPa) 

37.2024 

6.112 

36.91 

49.6032 

7.910 

43.88 


For 90*^ specimens, the value of SIF could not be evaluated. 
The reason for this was that the low strength of composite in 
transverse direction. Hence, the specimens fractured at loads 
which were so low that the maximum fringe order hardly exceeded 

oiae. 

Using these values of SIF and Eqs. (4.4). the stress field 
near the crack-tip was determined. The representative plots of 
OL. c^T and xlt with respect to the distance from the crack-tip 
have been shown in Figs. 4.5^, 4.?”, 4.7- 



45 - 





P'ot of longitudinal stress v/s 
distance from crock tio. 


£1 



ce From Crack (mm) 

Fig. 4.6 
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5 CONCLUSIONS: 

Increasing importance of fiber composites in various 
technological applications, calls for an accurate and efficient 
fracture analysis. For this purpose, the method of transmission 
photoelasticity can be an extremely powerful and reliable tool. 
In the present thesis, principles of photo -orthotropic elasticity 
have been employed successfully to understand the phenomena of 
fracture in composites. In this context, the following 
conclusions are in order: 

I- -A number of formulations for the fabrication of 

transparent composites as model materials for 

photoelastic study of composites have been reported in 
the literature. However, their properties especially 
transparency were not sufficient for the present 

investigations. The need of such materials is all the 
more acute if study of dynamic fracture in composites 
is intended. In the present investigations, a new model 
material for photo-orthotropic elastic studies has been 

developed. The optical properties of this model 

material as compared to those of other model materials 
are much better. The properties of the model are: 
Transmission ratio = 69% 

FcyL = 143.67 N/mm 
Fg’T = 66.67 N/mm 
FcrLT “ 54,07 N/mm 

2. The properties mentioned above can be modified by 

changing the volume fraction of the composite. 

Using this model material, and by employing principles 


3. 
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of orthotropic photoelasticity and orthotropic linear 
elastic fracture mechanics, the value of Stress 
Intensity Factor at different loads has been determined 
successfully. 

4. Comparing these values of S.I.F with the initial guess, 
it was inferred that the the values of stress intensity 
factor calculated as explained earlier, were reasonably 
correct . 

5. Using this value of S.I.F and making use of stress- 
field equations, the variation of cr^’ ^y '^xy with 
respect to the distance from crack tip was plotted. 

6. Examination of these plots revealed the presence of a 
singularity near the crack tip. Further, it was found 
out, that reversal in the sign of shear stress occurs 
near the crack tip, when 8 increases from 0*^ to 90'^. 
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APPENDIX-I 


Listing of FORTRAN Program for Finding Kj For Composites: 


slff v|^ ^ U/ ^ ^ slf* \M/ ^ \1/ vlv 1^ ^ \|/ ^ vi^ vJ> 

(tJ^ /f\ Jfj^ /J% ^1% ,/^ /f \ ./yS ,/|X yf% /^ 4ly\ i?|fv /f> 

EVALUATION OF STRESS INTENSITY FACTOR 

^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^1/ S^ *4f>^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^4^ >4^ *j^*4f* ^ 

4^JS rfs /fSi jy^ /|^ ^pfc /f\ ^j\ ^ys yf\ /pi ,/fs t^K jy^ /jk yjs. ^ih .ii^ /p^ ^w ^yx >|s /f\ /y% /f^ >j\ /f\ i/pi 

IMPLICIT DOUBLE PRECISION {A-H,0-Z) 

DIMENSION DEL(41,4),THT(41),X(41),Y(41),GK(41),RAD(41),PAR(4) 
DIMENSION N(41),SQR(5,5),COL(5,5).TRANS(4.41),GKMIN(41,4),SOL(4) 

C DEFINE THE FUNCTIONS 

PI=3.2429 
FL=i43.67E03 
FT=66.59E03 
FLT=54.067E03 
RNR=0.1000 
RTHTT=-0.001 

OPEN (UNIT=51,FILE=’FINAL.IN’ ) 

OPEN (UNIT=21,FILE=’SIFF.0UT’ ) 

READ(51,*)IA,P.B 

READ ( 5 1 , =♦= ) NPT . T , X0 . Y0 , REALKl , SIG0X 

STRESS=P/ ( B>«^T ) /1 . 0E+06 
WRI TE { 6, 1)IA, STRESS 
WITE( 21, 1)IA, STRESS 

1 FORMAT! ’Crack length : ’,12, ’ mm.’,/, 

$ ’Applied stress (e) : ‘,F7.3) 

M=4 

ITR=0 

PAR(1)=REALK1 

PAR(2)=RNR 

PAR(3)=RTHTT 

PAR(4)=SIG0X 

V?RITE( 6 , 2 ) PARC 1 ) /1 . 0E06 , PARC 2 ) , RTHTT , C SIG0X/1 . 0E+06 ) 

WRI TE C 2 1 , 2 ) PAR C 1 ) / 1 . 0E06 , PAR C 2 ) , RTHTT , C S IG0X/ 1 . 0E+06 ) 

2 FORMAT! ’ Initial Guess for K_1 = ’,F6.4,’ MPa4m’,/, 

$ ’N_r = ’ ,F6.4,/, ’e_r = ’ ,F6.4,/, ’ 0 _ox = ’.F7.3,’ MPa’) 

DO 10 1=1, NPT 
READC51,>K)XCI),YCI),NCI) 

CALL RADIUSCI,X,Y,X0,Y0,RADD,THTT) 

RADCI)=RADD 

THTCI)=THTT 

10 CONTINUE 

11 RNR=PARC2) 

RTHTT=PARC3) 

SIG0X=PARC4) 

DO 20 1=1, NPT 
RADD=RADCI) 

THTT=THTCI) 
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GALL FUNG ( I , N , PAE, THTT , RADD , FL , FT , FLT , T , GK) 

DO 15 

GALL DELTCI, J,N, PAR, X,Y,FL, FT, FLT,T,GK, DELTA) 

DEL(I, J)=DELTA 
15 CONTINUE 

20 CONTINUE 

C GENERATE THE TRANSPOSE OF DEL MATRIX 

DO 25 1=1, NPT 
DO 30 J=1.M 
TRANS(J, I)=DEL(I.J) 

30 CONTINUE 

25 CONTINUE 

C MULTIPLY TRANS & DEL TO GET 4 X 4 MATRIX 

CALL MATMPY (TRANS, DEL, SQR.M,NPT,M) 

C EVALUATION OF [ TRANSPOSE ^ COLUMN VECTOR 3 

DO 40 1=1, NPT 
GKMINCI, 1)=-GK(I) 

40 CONTINUE 

CALL MATMPY < TRANS, GKMIN, COL, M,NPT,1) 

C SOLVING FOR CORRECTIONS IN THE VALUES OF Ki etc. 

DO 43 1=1, M 
SQR(I,M+l)=COL(I,l) 

43 CONTINUE 

MP1=M+1 

CALL CHLSKY(SOR.M,MPl,SOL) 

RMAX=ABS ( SOL ( 1 ) /PAR( 1 ) ) 

DO 50 1=1, M 
PAR{I)=PAR(I)+SOL(I) 

RATIO=ABS ( SOL ( I ) /PAR( I ) ) 

IF (RATIO .GE. RMAX) RMAX=RATIO 
CONTINUE 
ITR=ITR+1 
WRITE(6,59)ITR 
WRITE(21,59)ITR 

WRITEC 6 , 60)PAR( 1 ) /1 . 0E06 , PARC 2 ) , PAR( 3 ) , PAR( 4 ) /1 . 0E+06 
WRITE(21,60)PAR(1)/1.0E06,PAR(2) ,PAR(3) ,PAR(4)/1 . 0E+06 
IF ( RMAX .LE. 1.0E-04 ) GOTO 55 
GOTO 11 

WRITE{ 6 , 60)PAR( 1 ) /I - 0E06 , PARC 2 ) , PARC 3 ) , PARC 4 ) /1 . 0E+06 
WRITE ( 21 , 60)PAR( 1) /1 . 0E06 , PARC 2 ) , PARC 3 ) , PARC 4 )/l . 0E+06 
FORMAT C ’ Iteration no . : ’ , 1 2 , / ) 

FORMATC/, ’K_l = ’,F14.6,’ MPa4ia’ , /, ’N_r = ’,F6.4, 

$ /,’S_r = ’ ,F6.4,/, ’a_ox = ’,F7.3,’ MPa’) 

STOP 
END 
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>>;)>: >K )K >l< >tc 5^: )(c >l< ^ >!< ;jc ^ >K 5)« * in * >|i: ^ ^ )(c ;K He ^ ^ ;j<;X ^ ^ jj; ^ ^ 

SUBROUTINE FOR THE VALUE OF RAD & THETA 

SUBROUTINE RADIUS ( I . X, Y , X0, Y0, RADD , THTT ) 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

DIMENSION X{41),Y{41) 

PI=3. 1429 


RADD-SQRT( (X{ I ) -X0) { Y( I ) -Y0) W ) 

IF ( (X(I)-X0) .EQ. 0) THTT=0.0 
IF ( (X(I)-X0) .EQ. 0) GOTO 1 
TI1TT=ABS(ATAN((Y{I)-Y0)/{X(I)-X0))) 

IF ({X(I) .LT. 0) .AND. (Y(I) .GE. 0)) 

IF ((Xd) .LT. 0) .AND. (Yd) .GE. 0) ) 

IF ((Xd) .GT. 0) .AND. (Yd) .LT. 0)) 

IF ((Xd) .GT. 0) .AND. {Y(I) .LT. 0)) 

IF ((Xd) .LT. 0) .AND. (Yd) .LT. 0)) 

CONTINUE 
RETURN 
END 


THTT=PI/2.0+THTT 
GOTO 1 

THTT=0.0-THTT 
GOTO 1 

THTT=PI+THTT 


SUBROUTINE FOR THE VALUE OF Gk 

SUBROUT INE FDNC (I , NN , XX , THTT , R , EL , ET , ELT , FF , G ) 
IMPLICIT DOUBLE PRECISION (A-H.O-Z) 

DIMENSION NN(41) .G(41),XX(4) 

PI=3.1429 

B2=l. 402818 

B1 =0.3056674 

F=1 . 0/SQRT ( 2 . 0>KPI!«R) 

C=Bl5t'B2*F/(B2-Bl) 

D=F/(B2-B1) 

ROl = ( COS ( THTT ) ^5jc2+Bl * ( ( SIN ( THTT ) ) . 25 

R02=( COS ( THTT )>t^^2+B2^{ (SIN (THTT) )^»^2) )*5^:0. 25 
GAl=ATAN(Bl>i^TAN(THTT) ) 

GA2=ATAN(B2*TAN(THTT) ) 

SIl=SIN(GAl/2.0) 

SI2=SIN(GA2/2,0) 

COl=GOS(GAl/2.0) 

CO2=COS(GA2/2.0) 

SIGMAX=XX( 1 ) W ( B2*C02/(R02 ) -B1*C01/(R01 ) ) +XX( 4 ) 

SIGMAY=XX(1)*D*(B2W1/(R01)-B1*G02/(R02)) 

T0WXY=XX(1)>KC^(SI1/(R01)-SI2/(R02)) 

AA=XX( 2 ) WS( 2 . 3 ) ) /FF 

BB=XX ( 2 ) =frSIN ( 2 . 0^XX( 3 ) ) /FF 
A= ( SIGMAY/EL-SIGMAX/ET+AA)=»'»<2 
B= ( 2>)fTOWXY/ELT+BB)=<'^2 
G (I ) =SQRT ( A+B) -NN (I ) /FF 
RETURN 
END 
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C SUBROUTINE FOR THE VALUE OF DIFFERENTIALS 

/i\ /js /|%, /p. \ yfv >p\ yj\ ^ ^ ^ \ \ ^ /f% ^ ^ ^ ^ jft. ^ yjs ^ ^ ^ ^ /f<i, ^ ^ ^ ^ ^ 

SUBROUTINE BELT ( I , J , NN , PAR , X , Y , FL , FT , FLT , T , GK , DELTA ) 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

DIMENSION NN(41),PAR(4),X(41),y(41),GK(41),GKK(41),PR(4) 
MM=4 

DO 10 ITR=:1,MM 
PR(ITR)=PAR(ITR) 

10 CONTINUE 

PR{J)=(1-1.0E-06)*PR{J) 

X00=0. 454545476135205D-003 
Y00=0. 454545476135205D-003 
CALL RADIUS { I , X , Y , X00 . Y00 , RADD , THTT ) 

CALL FUNC( I , NN , PR, THTT , RADD , FL , FT , FLT , T , GKK) 

DELTA= ( GK ( I ) -GKK ( I ) ) / ( 1 . 0E -06>i=PAR ( J ) ) 

RETURN 

END 

€ SUBROUTINE FOR MULTIPLICATION OF MATRICES 

SUBROUTINE MATMPY ( A1 . B1 . Cl , M , NNN , L ) 

IMPLICIT DOUBLE PRECISION (A-H.O-Z) 

DIMENSION Al(4.41) .B1(41,4),C1{5,5) 

DO 2 1=1, M 

DO 2 J=1,L 

C1(I,J)=0.0 
DO 1 K=1,NNN 

C1(I,J)=C1{I,J)+A1(I,K)!<'B1(K,J) 

1 CONTINUE 

2 CONTINUE 
RETURN 
END 

C SUBROUTINE CHLSKY 

SUBROUT INE CHLSKY ( A , N , M , X ) 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

DIMENSION A(5,5).X(4) 

C CALCULATE THE FIRST RON OF UPPER UNIT TRI. MATRIX 

DO 3 J=2,M 

3 A(1,J)=A(1,J)/A(1,1) 

C EVALUATION OF OTHER ELEMENTS OF U & L MATRICES 

DO 8 1=2, N 
J=1 

DO 6 II=J,N 
SUM=0.0D00 
JM1=J-1 
DO 4 K=1,JM1 

4 SUM=SUM+A{II,K)>»^A(K,J) 

5 A(II,J)=A(II.J)-SUM 
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IP1=I+1 
DO 7 

SUM=0.0D00 

DO 6 K=1,IM1 
SUM=SUM+A (I,K)*A(K,JJ) 

A(I, JJ)=(A(I, JJ)-SUM)/A(I, I) 
CONTINUE 

C SOLVE FOR X(I) BY BACK SUBSTUTION 

X(N)=A(N.N+1) 

L=N-1 

DO 10 NN=1,L 
SUM=0.0 
I=N-NN 
IP1=I+1 
DO 9 J=IP1,N 

3 SUM=SUM+A(I,J)>»^X(J) 

10 X(I)=A(I,M)-SUM 

RETURN 
END 



